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Abstract 

We study the topology dependence of the finite size corrections to the 
Ising model partition function by considering the model on a triangular lat- 
tice embedded on a genus two surface. At criticality we observe a universal 
shape dependent correction, expressible in terms of Riemann theta functions, 
that reproduces the modular invariant partition function of the corresponding 
conformal field theory. The period matrix characterizing the moduli param- 
eters of the limiting Riemann surface is obtained by a numerical study of the 
lattice continuum limit. The same results are reproduced using a discrete 
holomorphic structure. 
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I. INTRODUCTION 



In the context of two dimensional critical phenomena, finite size corrections 1 provide a 
direct connection between the critical behavior of specific lattice models and the conformal 
field theory description of the corresponding universality class. For a two dimensional system 
without boundaries, of typical length L, the partition function is expected from very general 
arguments 2 to behave, at the critical point and in the large L limit, as 

Z(L) ~ A L h exp (-/ooL 2 ) (1.1) 

where h depends on the system topology and singularities of the metric and A depends on 
the shape of the system. 

In this paper we study the shape dependence of the constant term A for the Ising 
model, and clarify its connection with conformal field theory, by considering the model on 
triangular lattices embedded on a genus two surface. Finite size corrections on higher genus 
lattices have been previously studied on square lattices 3 for both the Ising model and the 
close-packed dimers. Here we test the universality of the results obtained by considering 
the triangular lattice case. The dimer model is not considered since it is not critical on 
triangular lattices 4 . 

For a lattice embedded on a genus g surface, the Ising model partition function can be 
expressed, using the Kasteleyn dimer formalism 5-9 , in terms of Pfaffians of 4 9 adjacency 
matrices A k 

Z g (N i ,K i ) = '£PfA k (N i ,K i ). (1.2) 

k=i 

where the iVj and Ki are respectively the integer sizes and coupling constants characterizing 
the lattice. We study numerically these Pfaffians in a number of genus g — 2 triangular 
lattices and show that they satisfy at the critical point and in the large lattice size N limit 



PfA k (N t ,K t ) 
A/-™oc Pf MN^KA 



e[k\(o\n) 



;i. 3 ) 



where the Q[i](0|fi) are genus two Riemann theta functions, defined in Section III, and Q is 
the 2x2 period matrix. The shape dependent term in the partition function (1.1) is found 
to be, in the large lattice limit, 



16 



A^A £|e[*](0|n)| (1.4) 



k=i 



where A is a /c-independent constant. 

This sum over theta functions can be identified with the classical winding part of the 
partition function of the corresponding conformal field theory 10-12 on a genus two Riemann 
surface with a complex structure characterized by the period matrix Q. 

The remarkable similarity between the combinatorial result (1.2) and the sum over A 9 
fermion boundary conditions 10 in conformal field theory, together with the explicit results 
obtained so far for genus one 13 ' 14 and two 3 suggest that the shape dependent correction (1.4) 
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is universal. By this we mean the following; the expression (1.4), with a sum over 4 9 theta 
function for more general topology, is found at criticality for all lattices independently of the 
lattice structural details, such as aspect ratios and coupling constants, while these structural 
details determine the g x g period matrix Q. 

We are left with the problem of understanding the dependence of the period matrix 
Q on the lattice detailed structure. For toroidal lattices this is simply the lattice aspect 
ratio weighted by the coupling constants along the two lattice directions 13-16 . In the higher 
genus case we will show that the lattice period matrix Q can be evaluated as the continuum 
limit of a discrete period matrix obtained by solving certain sets of finite difference Poisson 
equations on the lattice. 

The concept of discrete period matrix was first introduced by the authors 3 , in the con- 
text of locally square lattices, using a discrete formulation of harmonic and holomorphic 
differentials. While discrete harmonic differentials are well understood, in the context of 
combinatorial Hodge theories on lattice cochains 17-19 , discrete holomorphy is more elusive, 
the main obstacle to a lattice formulation being the discretization of the Hodge operator. 
Discrete Hodge operators have a long history in finite difference calculations in electrody- 
namics and other problems of computational physics 20-25 . Two dimensional formulations 
have also been studied 26 ' 27 in connection with the Ising criticality and discrete holomorphy. 

In this paper we extend our previous study 3 to triangular lattices. We show that the 
period matrix for a triangular lattice G can be determined in terms of an effective square 
lattice Gn- This correspondence depends explicitly on the Ising model criticality condition 
and provides a geometrical interpretation of this condition. 

The structure of the paper is as follows: 

In section II we discuss the Kasteleyn formalism and define the adjacency matrices in 
(1.2) for the genus two lattice shown in Fig. I. In section III we study (1.3) numerically 
by evaluating ratios of determinants of adjacency matrices. In section IV we define discrete 
harmonic differentials. In section V we define discrete holomorphic differentials and the 
lattice period matrix and discuss the relation between our formalism and the formalism 
of reference 27 . The period matrices of (1.4) are reproduced using discrete holomorphy in 
section VI. Finally, our conclusions are presented in section VII. 

II. THE KASTELEYN FORMALISM 

The partition function of the Ising model on a lattice, or more generally on a graph, em- 
bedded without superposition of edges on an orientable surface of genus g can be expressed as 
a sum (1.2) over the Pfaffians of 4 9 adjacency matrices. This result, due to Kasteleyn 7 ' 28 ' 8,9 , 
is the starting point for the study of the Ising model on higher genus lattices. 

In this section we present the main aspects of the Kasteleyn formalism and define the 
adjacency matrices for the genus two triangular lattice, G, shown in Fig. 1. This lat- 
tice is characterized by five integers sizes (Mi, M 2 , K, Ni, N 2 ) and three coupling constants 
(Kh,K v ,Kd) distinguishing the horizontal, vertical and diagonal edges. The boundary condi- 
tions are such that the lattice can be drawn without superposition of edges only on a surface 
of genus two or higher. The cycles of lattice edges a i: b i: with % — 1, 2, represent a canonical 
basis of the first homology group of the embedding surface. The same cycles drawn over the 
dual lattice edges will be denoted by di, bi (see Fig. 2). 
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The Ising model on the lattice G can be expressed in terms of a dimer model on a 
decorated lattice G' using a well known 6 ' 29 one to one correspondence between polygon 
configurations of the Ising model high temperature expansion in G and close-packed dimer 
configurations in G' . A close-packed dimer configuration is a selection of lattice edges such 
that every vertex is included once and only once as a boundary of a selected edge. 

The decorated lattice G' is obtained by replacing each vertex of the original lattice G 
with a decoration graph. For a triangular lattice the decoration can be chosen 29 to be a 
hexagon with all vertices connected (see Fig. 2) and the decorated lattice will have 6J\f 
vertices for a original lattice of size Af. 

Each edge on the decorated lattice G' is assigned a weight: Wi = tanh f3Ki for edges of 
type % inherited from the lattice G and a weight 1 for the internal edges of the decoration. 
The Ising partition function on the lattice G can then be expressed as 

G' 

Z Mng = (2 cosh K v cosh K h cosh K d ) N £ w v nv w h nh w d nd (2.1) 

dimer config. 

were the sum runs over all the close-packed dimer configurations of the decorated lattice G' 
and rii is the number of edges of type i in a given dimer configuration. 

The dimer partition function can be expressed in terms of Pfaffians of the Kasteleyn 
adjacency matrices. These are signed adjacency matrices defined in the following way: 
label the decorated lattice vertices with an integer from 1 to M. and choose a lattice edge 
orientation by assigning to each edge a direction represented by an arrow (see Fig. 2). The 
signed adjacency matrix corresponding to this edge orientation is the M. x M. matrix 
with entries 

!z if there is an arrow from vertex i to vertex j of weight z 
—z if there is an arrow from vertex j to vertex i of weight z (2.2) 
otherwise 

The Pfafnan of such an M. x M anti-symmetric matrix, with M. even, is defined as 

Pf (A) = 2_m/2^^^/2)! ^ €p ^pi-pz^pzP' 1 ' ' ' Apm-ipm (2-3) 

where the sum goes over all the permutations p of the integers from 1 to M. and e p — ±1 
for even and odd permutations respectively. 

From the definition of the adjacency matrix it is clear that each non zero term in the 
Pfafnan equals, in absolute value, a term in the dimer partition function, with the relative 
sign depending on the choice of the edge orientation. If we choose an orientation such that 
all the terms in the Pfafnan have the same relative sign, the Pfafnan of the adjacency matrix 
equals the partition function modulo an overall sign. 

Kasteleyn showed in a beautiful tour de force of combinatorics 5 that edge orientations 
with this property exist. These are the edge orientations such that every lattice face has an 
odd number of clockwise oriented edges. In a genus g lattice there will be 4 9 inequivalent 
such orientations and the dimer partition function is given as a linear combination of their 
Pfaffians, see references 7-9 for more details. 

For the genus two decorated lattice G 1 the relevant sixteen edge orientation classes can be 
labelled as A(n 5l , n^, n 52 , n~ b2 ) with n x — 0, 1 for x = 01,02,61,62- The orientation A(0000) 
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is shown in Fig. 2. An orientation with a certain n x = 1 is obtained from the corresponding 
orientation with n x = by introducing a disorder loop along the non-trivial cycle x, this is, 
by reversing the orientation of all the edges crossed by the cycle x. 

To make connection with the theta functions characteristic and allow for more compact 
equations we will also use the alternative notations 

= Ai (2.4) 

with the integer label given by % = 16 — 8n 5l — 4n^ — 2n 52 — n~ h ^. In terms of these adja- 
cency matrices the Ising model partition function on G is given by 

Z = \{ P 1 - P 2 - P 3 - P 4 - P 5 + P 6 + P 7 + P 8 - 
P 9 + P w + P u + 

Pi 2 — Pl3 + P 14 + Pl 5 + Pw ) ' ' 1 

where Pj = Pf Ai and the adjacency matrices are defined on the lattice G'. 

For translational invariant lattices, such as the torus triangular lattice, the Pfaffians of 
the adjacency matrices can be evaluated in a closed form and the theta function dependence 
extracted by a careful asymptotic analysis 13,14 For a genus two lattice as G' such an analytic 
treatment is not possible and we are forced to resort to numerical evaluations of the Pfaffians. 



^.(n^, n^, n fl2 , n b2 ) — A 



III. RATIOS OF DETERMINANTS AND THETA FUNCTIONS 



The Pfaffians of the adjacency matrices defined on the previous section can be nu- 
merically evaluated using the fact that Pf A = Vdet A for an anti-symmetric matrix A. 
We evaluate these determinant for a sequence of lattices characterized by the integers 
sizes (Mi, M2, K, Ni, N%) = (mi, 777,2, k, ni, 712) L with increasing L and fixed aspect ratios 
(mi,m2, k, Tii, 712). The Ising model coupling constants were restricted to the ferromagnetic 
case Ki > and chosen to satisfy the criticality condition 



smh(2K h ) sinh(2X„) + sinh(2iY fe ) sinh(2iY d ) + sinh(2ir„) sinh(2iY d ) 



1. 



(3.1) 



In order to eliminate the common bulk and algebraic terms, we consider ratios of the 
sixteen determinants to the largest among them 



Ri — 



det Ai 
det A 

TT1>V 



(3.2) 



An example of the results obtained is shown in Fig. 3 and Fig. 4, where the fifteen ratios 
are shown, for a lattice with aspect ratios (777,1, 777,2, 77,1, 77,2) = (4,2,2,2,4) and coupling 
constants [wh, w v , Wd] = [0.369, 0.159, 0.282]. The ratios are shown as function of the number 
of lattice vertices, M . Six of them, shown in Fig. 4, are found to vanish in the large M limit 
while the remaining nine converge to finite values, Fig. 3. 

We found that these ratios in the large M limit can be expressed in terms of ratios of 
genus two Riemann theta functions. These functions are defined 30 by the quickly converging 
series 



e 



OL 



(z,fi) 



neZ' 2 



i7r(n + ct) T Q(n + a) + 2vn(n + a) T (z + j3) 



(3.3) 
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where a, (3, z and n are 2-vectors half-integers, complex numbers and integers respectively. 
The 2x2 period matrix Q is a symmetric complex matrix with positive definite imaginary 
part 



^11 ^12 
^12 ^22 



+ i 



^11 ^12 
^12 ^22 



(3.4) 



The extrapolated values of the determinants ratios, in the L — > oo limit, can be expressed 



as 



det [A 


Cl C 2 
C3 C4 


]) 


det (A 


'00" 



) 






' ci/2 c 2 /2 " 
. c 3 /2 c 4 /2 _ 


(0,Q) 


e 


"00" 
_ _ 


(0,Q) 



(3.5) 



for the 16 combinations of q = 0, 1, with a remarkable precision from to 10 -6 , with a 
period matrix Q determined by a suitable numerical fitting procedure. Examples for three 
different lattices are shown in table I. For each lattice the first column shows the L — > oo 
extrapolated ratios of determinants (3.2) and the second column shows the theta function 
ratios 



@(16-8di-4 Cl -2ci 2 -C2)(^) — © 



ci/2 c 2 /2 
dl/2 <2 2 /2 



(o,n)/e J J (o,n) 



(3.6) 



for Ci,di = 0,1, with the fitted period matrix. 

The period matrices in (3.5) are in general complex for triangular lattices unlike what 
was found for square lattices 3 where the period matrix is purely imaginary. Away from 
criticality, the sixteen determinants converge rapidly to the same value and (3.5) does not 
hold. 



IV. HARMONIC DIFFERENTIALS ON THE LATTICE 

We now turn to the problem of evaluating the period matrix Q in (3.5) directly from the 
lattice structure using a discrete formulation of harmonic and holomorphic differentials. 

In continuum Riemann surface theory the period matrix characterizing a given surface 
is evaluated in two steps, the first step being the determination of the space of harmonic 1- 
forms on the surface, the second step being the decomposition of this space into holomorphic 
and anti-holomorphic sub-spaces. On a lattice the role of p-forms is played by p-cochains, 
which are linear functionals defined on formal sums of the lattice p-elements: vertices, edges 
and faces. We will refer to the p-cochains by lattice functions, lattice differentials and lattice 
volume forms. 

In this section we give the first step towards a definition of the lattice period matrix 
by defining and evaluating the harmonic lattice differentials or 1-cochains on the lattice. 
Discrete harmonicity and the relation between p-cochains and p-forms is well understood 
from the work of Eckmann 17 and Dodziuk and Patodi 18,19 on combinatorial Hodge theories. 
Here we generalize the treatment done in reference 3 to triangular lattices by showing that 
harmonic differentials on the triangular lattice are completely determined by their restriction 
to an effective square lattice. 
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A. Lattice differentials on a triangular lattice 



A lattice function / is determined by its value on the lattice vertices f[n]:n = l,...,Af 
with n being an integer labeling of the M vertices in the lattice. 

A lattice differential w is determined by its value on the lattice edges. For the lattice G 
shown in Fig. 1 this dependence can be expressed as w[n\k] : n — 1, . . . ,Af; k = 1, 2, 3 where 
[n|l] stands for the horizontal edge immediately right of vertex n, [n\2] for the vertical edge 
immediately below and [n|3] for the diagonal edge laying in between, see Fig. 5. Edges are 
oriented as shown in the figure according to their classification into horizontal, vertical or 
diagonal. The integral of a lattice differential w along a path C of lattice edges is defined 
to be the sum of the values of w on all the edges e included in the path 



w 



£±Me] (4.1) 



C e&C 

the minus sign applying to edges with opposite orientation to the one of the path. 

A lattice volume form 77 is determined by its value on each lattice face. The genus 
two triangular lattice G has 2J\f — 2 faces, two of which are hexagons the remaining being 
triangles. A volume form 77 in G can be expressed as rj[q] : q = 1, . . . , 2J\f — 2 where q is an 
integer labelling of the faces of G. The integral of a volume form 77 over a given area A on 
the lattice is the sum of the values that 77 takes on all the lattice faces q contained in that 
area 



»7=5>[g]. (4.2) 



The lattice exterior derivative d is a linear operator that takes lattice functions into 
lattice differentials and lattice differentials into lattice volume forms. It can be defined as 

(df)[e] = f[n el ]-f[n e2 ] (4.3) 
(dw) [g] = £±Me] 

egg 

where n e \ and n e 2 are the two vertices corresponding to the end and beginning of the 
oriented edge e and the sum runs over all the edges e in the boundary of the face q, the plus 
or minus sign being determined by the relative orientation of the edge e to the anticlockwise 
orientation of the face. 

The exterior derivative defined in this way satisfies a discrete version of Stokes theorem. 
Let C (n, n') be a path of lattice edges from vertex n to vertex n' and A an area on the 
lattice, then we have that 

d f = f[n'] - f[n] (4.4) 

C(n,n') 

dw = w 

A JdA 

where OA is the path along the boundary of the area A with an anticlockwise orientation. 
An inner product on lattice functions, differentials and forms can be defined as 
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(/,/') = 5>(n) /[n]/'[n] (4.5) 



(tu, u/) = Pi(e) tw[e] iu'[e] 



(v,v') = J2p^) v[q]v'[q] 

q 

where the pi are positive functions on the vertices, edges and faces. For our purposes, in 
the lattice G we define po(n) = pi{q) = 1 and let pi(e) take one of three values pn, Pv 
and pd for horizontal, vertical or diagonal edges respectively. In terms of the edge labeling 
introduced above we have then 

(w, w') = ^2 (j>h w[n\l] w'[n\l] + py w[n\2] w'[n\2] + po w\n\?>\ w'\n\S^ (4.6) 

n 

From a choice of inner product there follows a definition of the lattice co-derivative 5, 
the adjoint operator of the exterior derivative, 

(w,df) = (5wj) (4.7) 
(rj,dw) = (5r],w). 

The co-derivative takes lattice differentials into lattice functions and lattice volume forms 
into lattice differentials and can be expressed as 

(5w)[n] = J2±Pi(e)w[e} (4.8) 

Pl\ e ) 

(4.9) 

where the sum on the first equation runs over all edges e that contain the vertex n, the plus 
and minus sign corresponding to n being the ending or the starting vertex of the oriented 
edge. On the second equation q d and q er stand for the faces left and right of the edge e with 
respect to its orientation. For the lattice G we have explicitly 

Sw[n] = p H (w[left(n)|l] - w[n\l\) + p v (w[n\2] - w[up(n)\2}) (4.10) 
+p D (w[\eh(n)\3] -w[up(n)|3]) 

where the functions right(n), left(n), up(n) and down(n) give the label of the vertex immedi- 
ately right, left, above and below of vertex n, see Fig. 5. 

Given a definition of the lattice exterior derivative and co-derivative we define, in exact 
analogy with the continuum concepts, a lattice differential w to be closed if dw — 0, exact 
ilw — df, co-closed if 5w = and co-exact if w = Srj. A lattice differential is said to be 
harmonic if it is both closed and co-closed 

. , . | dw\q] = 0, for all faces q in G , . , lN 

w is harmonic = < r r L i _ „ „ , . . „ (4.11) 
ydw[n\ = 0, tor all vertices n m G. 

These equations can be seen as linear system of equations on the unknowns We 
will see that for a genus g lattice there are 2g linear independent solutions, corresponding 
to the 2g dimensional vector space of harmonic differentials of the continuum theory. 
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B. The effective square lattice description 



The linear system (4.11) can be partially solved by noting that the value of a closed 
differential on the diagonal edges w[n|3] is determined by its values on the adjacent horizontal 
w[n|l] and vertical tu[n|2] edges. Using this constrain to eliminate the w[n|3] from the 
second equation in (4.11) the problem of determining the lattice harmonic differentials on 
the triangular lattice becomes defined on an effective square lattice. 

Let Gn be the square lattice obtained from the triangular lattice G by removing all the 
diagonal edges. Consider Fig. 5; any square face q of G u contains two triangular faces of G, 
t qi and t q2 . These two triangular faces share the diagonal edge [n|3] that enters the closeness 
condition in (4.11) in two equations 

dw[t qi ]=0, dw[t q2 } = (4.12) 

that can be equivalently expressed as 

d n w[q] = (4.13) 
w[n\3] = w[n\l] + w[n\2] (4.14) 

where d a is the exterior derivative on the square lattice G a and the second equation deter- 
mines the value of the differential on the diagonal edge. 

In a similar way the hexagonal faces of G can be seen as contained on an octagonal face 
of Gn together with two triangular faces, see Fig. 5, The closeness condition on these three 
faces of G is 



dw[t 1 } = 0, dw[t 2 ] = 0, dw[t 3 } = (4.15) 

and can also be expressed in terms of the vanishing of exterior derivative do on the octagon 
q and condition (4.14) for the two diagonal edges involved 

d D w[q] = (4.16) 
iu[n|3] = + iu[n|2] for n — m,m'. 

Using (4.14) to eliminate the tu[n|3] terms in the co-closeness condition in (4.11) we 
obtain that the harmonic differentials on the triangular lattice G are completely determined 
by the following conditions on their restriction to the squared lattice G u 

. , . [ d a w\q] — 0, for all faces q in G n , A ,„ N 

w is harmonic = < r r „ , L . . „ 4.17 

I o n w[n\ = 0, tor all vertices n m G D 

where the finite difference operators on Gn are given by 

daw[q] = £(w[down(&)|l] + ^[right(^)|2] - w[qi\l] - w[qi\2]) (4.18) 

i 

5 a w [n] = p H (w[left(n)|l] - w[n\l\) + p v (w[n\2] - w[up(n)|2]) (4.19) 
+p D (iu[left(n)|l] + w[left(n)|2] - iu[up(n)|l] - w[up{n)\2]) 
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the sum in (4.18) goes over i — 1,2 for the two octagonal faces of Gn, with qi and q<i being 
the labels of the two vertices that can be seen as the octagon upper left vertices in Fig. 1; 
and over i — 1 for the squared faces of G n of which q\ is the upper left vertex. 

This notation 3 for the faces and vertices of a locally square lattice, has the advantage of 
relating the labeling of vertices with the labeling of faces: if q labels a face then label the 
upper left vertices of that face. Reciprocally if n labels a lattice vertex we define n to be the 
label of the face of which n can be seen as the upper left vertex. We then have the relations 

& = q (4.20) 
n E {h 1: h 2 } (4.21) 



C. Evaluation of the harmonic differentials 

For a given lattice G a basis of the space of harmonic differentials can be evaluated by 
direct solution of one of the equivalent linear systems (4.11) or (4.17). 

The number of solutions can be bounded from below by subtracting the number of 
unknowns by the number of equations. Two additional solutions are provided by the two 
dimensional space of constant differentials and the total number of solutions is at least 

#(edges) - #(vertices) - #(faces) + 2 = ~x + 2 = 2g (4.22) 

where x is the Euler characteristic of the lattice. There are no additional solutions because 
harmonic differentials are completely determined by its periods, or integrals, along the non- 
trivial homology cycles of the lattice and 2g is precisely the dimension of the first homology 
group (see 3 for further details). 

The linear system (4.17) can always be evaluated numerically, although the evaluation 
of the kernel of a large linear system is computationally demanding. For practical purposes 
there is a much more efficient method based on a discrete version of the Hodge decomposition 
theorem. This theorem states that the space of lattice differentials in G a has an orthogonal 
decomposition in terms of exact, co-exact and harmonic differentials and that any lattice 
differential w can be written in an unique way as 

w = d n f + h + 5oV (4.23) 

where h is a harmonic. The orthogonality of the different kinds of differentials follows 
directly from (4.17) and the property d n d n = 5a8a = 0. 

This result allows the determination of harmonic differentials by orthogonal projections. 
Our objective is to obtain a basis {A k , B k : k — 1, 2} of the space of harmonic differentials 
satisfying the normalization conditions 

(4.24) 



/ A k 


— S k j, 


[ A k = 






Jb 3 


I B k 


= o, 


l B k = 5 k j 
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with the dj, bj being any choice of closed paths on the lattice representing a basis of the first 
homology group. We proceed in the following way: start with closed but not harmonic dif- 
ferentials A k , B k with the periods required by (4.24). A possible choice for these differentials 
is shown in Fig. 6 where we take A k (B k ) to be zero on all edges except the ones crossed 
by the dual lattice cycles b k (respectively a k ). Since closed differentials are orthogonal to 
co-exact differentials it follows from (4.23) that closed differentials can be written as the 
sum of a harmonic differential with the same periods and an exact differential 

A k = A k + d u f£ (4.25) 
B k = B k + d n f^ 

applying 5 on the these equations and using (4.17) we obtain 

/\nf k A = 5nA k (4.26) 
A n f* = 6 n B k 

where the Laplacian operator A n = 5 a da on lattice functions is given explicitly by 

(A n f)[n] = 2{p H + Pv + p D ) f[n] (4.27) 
-p H (/[right(n)] + /[left(n)]) 
~Pv (/[up(n)] + /[down(n)]) 
-p D (/[right(up(n))] + /[down(left(n))]). 

The Laplacian operator can be seen as a Af x Af matrix of rank (Af — 1) acting on the 
functions A/"-vector. The linear system of equations (4.26) can then be solved numerically for 
reasonably large Af by fixing the value of the functions at a lattice vertex and the solutions 
can be differentiated and subtracted from the A k , B k to obtain a normalized basis of the 
space of harmonic differentials {A k ,B k : k = 1,2}. This method is numerically more 
efficient than the direct solution of the linear system (4.17). 

V. HOLOMORPHIC DIFFERENTIALS AND THE PERIOD MATRIX 

The second step in the definition of a lattice period matrix is the decomposition of the 
space of lattice harmonic differentials, obtained in the previous section, into holomorphic 
and anti- holomorphic sub-spaces. For general values of the edge weights p#, py and Pd in 
(4.6) this is not possible. However if the edge weights are related through the Ising model 
criticality condition we are able to interpret dn and Sn as the geometrical realization of a 
square lattice tilted by an angle and a generalization of the square lattice formalism 3 allows 
the definition of a discrete period matrix for the triangular lattice G. 

A. Geometrical interpretation of the criticality condition 

In order to reproduce the period matrix in the Ising model finite size correction (1.4) the 
edge weights p#, py and p^ must be given in terms of the Ising model coupling constants 
along the same direction as 
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Pj = sinh2.Kj. 



(5.1) 



Since the coupling constants at criticality verify (3.1) we can parameterize the edge weights 
Pi in terms of a positive number r and an angle n/2 < 9 < n 

Ph = -^—(r + cosO) (5.2) 
smv 

1 ( l 
Pv = — — j. - + cos 9 



sin 9 \r 
— cos 9 

sin 9 

with 



- cos 9 = PD (5.3) 

\J(pH + Pd)(pv + Pd) 



r 



. pL±££>. (5.4) 
\ (Pv +Pn) 



This parameterization has a clear geometrical interpretation with r being the ratio of 
length scales along the vertical and horizontal directions and 9 a tilt angle, see Fig. 7. To 
show this, consider the co-derivative (4.19) in terms of the new parameters, dropping an 
overall sin 9 factor, we have that 

5 D f[n] = r (w[\eft(n)\l\ - w[n\l}) + 1/r (w[n\2] - w[up(n)\2\) (5.5) 
-cos6» (w[n\l] -w[up(n)\l] + w[\eh(n)\2] - w[n\2\) . 

and the discrete Laplacian operator on lattice functions (A n = 5 n d n ) is given by 

(A D f)[n] — 2(r + 1/r + cos 9) f[n] (5.6) 
-(r + cos^) (/[right(n)] + /[left(n)]) 
-(1/r + cos^) (/[up(n)] + /[down(n)]) 
+ cos0(/[right(up(n))] + /[down(left(n))]) 

Modulo a constant term, this Laplacian is the discretization of a continuum Laplacian 
describing a geometry with an angle 9 between the directions x and X\ 



with r being the ratio of the lattice spacings along these directions, precisely the geometry 
of Fig. 7. 

The Ising model criticality condition (3.1) allows the description of discrete harmonic 
differentials on a triangular lattice as harmonic differentials on a square lattice tilted by 
an angle 9. While the Ising criticality condition has been previously given a geometrical 
interpretation 26 at the level of discrete holomorphy, the reasoning above shows that a geo- 
metrical interpretation exists already at the level of discrete harmonicity. 
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B. The Hodge operator 




In the continuum theory the holomorphic and anti-holomorphic sub-spaces are obtained 
from the space of harmonic differentials by the projection operators 

P = (l + i*)/2 and P = (l-i*)/2 (5.8) 

where * is the (continuum) Hodge operator 31,32 . This operator is an endomorphism on the 
space of harmonic differentials and satisfies * 2 = — 1. Its action on differentials is given by 

* [f x (x, V) dx + f y (x, y)dy] = -f y (x, y) dx + f x (x, y)dy (5.9) 

and a differential is said to be holomorphic if it is of the form Pw with w harmonic. 

We will proceed in a similar way in the discrete theory and define lattice holomorphic 
differentials by a projection with a discrete Hodge operator *. This operator on lattice 
differentials should be defined in such way that (1 ± i*)/2 are projection operators and 
endomorphisms on the space of harmonic differentials, or equivalently that 

= 0, for all q j d * w[q] = 0, for all q , > 

= 0, for all n | 5 * w[n] = 0, for all n 

= -1 (5.11) 

There is no known definition of a discrete Hodge operator that satisfies these properties 
precisely while keeping the dimensions of the space of holomorphic differentials equal to g. 
We propose the following definition 

Definition: the Hodge operator on the effective square lattice G u with edge weights 
parameterized according to (5.2) is defined by 

1 1 cos 9 

(*w) Ml] = — - Mup(ra)|2] + — -w[up(n)\l] (5.12) 
r sin 9 sin 9 

(*w) [n\2]= r ^—w[\eh(n)\l] - ^ iu[left(ra)|2] 
sine' siny 

For 9 = 7i / 2 we recover the discrete Hodge star for the simple square lattice 3 , minor 
differences are due to the different labeling we have chosen for the lattice edges, see Fig. 9 
in that reference. 

This definition of the discrete Hodge operator relates the exterior derivative d with the 
co-derivative 8 in a way that closely satisfies (5.10). From the definition above, (4.18) and 
(5.5) it follows that 

(d™)[?] = --^£Mfc] (5-13) 
bin a i 

(6*w)M = \ " (n) if n e { n i> n 2, n 3 ,n 4 } ( , 

1 m \ sin0diu[up(left(n))] + ^{5w[up{n)\ - 5w[\eft(n)}) otherwise 1 ' 

where in the first equation the sum goes over i—1,2 for the octagonal faces and over i—1 for 
squared faces. In the second equation, the rik are the four lattice vertices where left(up(nfc)) ^ 
up(left(n fc )) and up(n) stands for the tilde of the vertex up(n). 
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Because of the four vertices the discrete Hodge operator fails to be an endomorphism 
on the harmonic lattice differentials: the requirement that a differential w is closed and 
co-closed is not sufficient to ensure that the corresponding *w is co-closed. This is due to 
the fact that the lattice has more vertices than faces and co-closeness is a condition verified 
on vertices. The condition is still partially satisfied in the sense that, if n\ and n 2 are on 
the same octagon, then the a{rik) satisfy 

ot{n\) + a(n 2 ) = sin 9 d iu[up(left(ni))] (5.15) 
cos 

+——{5w[up{n 1 )] - £iu[left(ni)] + 5w[up(n 2 )] - 5w[\eh{n 2 )}) 
sm.9 

and similarly for a(n 3 ) + 01(714). 

Unlike the continuum Hodge star, our lattice definition does not satisfy (5.11) but rather 

* 2 w[n\l] = -iu[left(up(ra))|l] + ^ ( D(*w)[n] + ^ D(w)[\eft(n)]\ (5.16) 

sm9 \ sin 6* J 

* 2 w[n\2] = -w[up(left(n))|2] + r—^-D{w)[\eh{n)] 

sin 9 

were D(w) is defined as 

D(w)[n\ = (iu[up(ra)|l] - w[\eft(n)\l]) - l - cos6 (w[up(n)\2] - w[\eft(n)\2}) . (5.17) 

as in the simple square lattice case 3 the action of * 2 on lattice differentials produces a 
diagonal translation, the size of a lattice spacing, that becomes negligible in the continuum 
limit. 



C. Differentials in the dual lattice 

We can also consider the construction of harmonic and holomorphic differentials on Gp, 
the dual lattice of Gn- While most of the concepts above follow through with minor changes 
there are a few relevant differences that make the dual lattice a testing ground for the 
dependence of our definitions on the lattice structure. 

We define G* n as the dual of G a in the usual graph theoretical sense together with the 
same assignment of edge weights, parameterized by r and 9, for horizontal and vertical edges. 
The two lattices Gn and G n differ at the conical singularities: where one has an octagonal 
face the other has a vertex with eight adjacent edges. This fact has important consequences 
in equations (5.13) and (5.14). As mentioned before closed differentials w in G a do not 
generate co-closed differentials of the form *w because the lattice has more vertices than 
faces. For the dual lattice Gp, where the number of faces is larger than the number of 
vertices, the opposite happens and (5.13) is given in terms of 0(71) terms not expressible in 
terms of derivatives and co-derivatives at a single point. 

For calculation purposes the dual lattice can be schematically represented as Gn with 
two vertices identified on each octagonal face , see Figs. 8 and 9. The faces of G* D are all 
squares and two of the lattice vertices have eight adjacent edges. Since the edges of G^ are 
in one to one correspondence with the edges of G u and we can use a edge labeling for G^ 
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similar to the one used for Gn- The Hodge operator and the exterior derivative on functions 
can be defined in the same way for both lattices while the exterior derivative d n and the 
co-derivative 5 D are defined as 

d* a w[q*\ = {w[q*\l]+w[r\ght{q*)\2\ - w[up{q*)\l] - w[q*\2\) (5.18) 
5* n w[n*]=J2PH (Mdown(0|l]-M<|l])+p y (w[n*\2] - w[up(n*)\2}) (5-19) 

i 

+ PD (w[\efr(n*)\l] + w[\efr(n*)\2] - w[n*\l] - w[n*\2}) 

where in the first equation q* is the vertex of G D that is the upper left corner of the square 
face q* of G* a . In the second equation the sum goes over % — 1, 2 when n* is one of the G D 
vertices with eight edges and i — 1 for all other vertices. The functions right(n), left(n), up(n) 
and down(n) and the relations between vertices, edges and faces in £?□ are evaluated on Gn 
with values in G^ which means that we start by distinguishing over different n* but identify 
them in the final result. This set of rules provides a complete prescription to evaluate 
holomorphic differentials in G a . 



D. The torus modular parameter 

As an illustration of the formalism developed above we will consider the genus one 
toroidal lattice and reproduce the modular parameter, known by exact asymptotic analysis 14 . 

Consider a toroidal triangular lattice with M rows and iV columns, characterized by 
the coupling constants Kh, K v and along the horizontal, vertical and diagonal edges. 
Let {a, b} be a basis of the first homology group with a being a horizontal loop and b a 
vertical loop. The harmonic differentials, solution of (4.17), on this lattice are the constant 
differentials 

A[n][l] = a, A[n][2] = (3 (5.20) 

for all vertices n, with a and f3 being two arbitrary constants. These harmonic differentials 
can be projected into the space of holomorphic differentials as 

(1 + i*) A[n] [l]= a + i(^?- a --J-p] (5.21) 

\sin0 r sin0 y 

(1 + i*) A[n\ [2] = (3 + i (-^ (5 + r J- a) (5.22) 

\ sin t) sin V J 

for all vertices n. The lattice modular parameter r is obtained by integrating this differential 
along a loop of the homology class of b after a proper normalization is chosen 

f(l+i*)A=l (5.23) 

J a 

f(l+i*)A = T (5.24) 

Jb 

immediately we reproduce the know result 14 

M 

r = — r (cos 9 + i sin 9) (5.25) 
with r and 9 given by (5.1) and (5.3). 
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E. The period matrix on higher genus lattices 

After obtaining a basis {A k , B k } for the space of harmonic differentials on a given lattice 
G, using the procedures of section IV, we proceed to evaluate a basis for the space of holo- 
morphic differentials of the form {T k = (1 + i*)H k : k — 1 . . . , g} where the H k are harmonic 
differentials, satisfying the normalization conditions 

f {H l + i*H l ) = 8 kl (5.26) 

Ja k 

[ (H l + t^H l ) = n kl . (5.27) 

Equation (5.26) is a set of 2g real constraints on H k from which we can determine the 
coefficients of H k in the basis {A k , B k } and equation (5.27) is our definition of the finite size 
lattice period matrix. 

The same construction can be made on the dual lattice G* a . We start by obtaining a 
basis of the harmonic differentials {A* k , B% : k — 1, 2} using the procedure of section IV and 
(5.18) and (5.19). Then using the same projection operator we obtain a basis for the space 
of differentials {Tl = (1 + i*)H k : k — 1 . . . , g} with the H k being harmonic differentials in 
Ctq, satisfying the normalization conditions 

/ (H; + i*Hl)=5 kl (5.28) 

Ja k 

[ (H; + i*H;) = n* kl (5.29) 

The discrete period matrix on the dual lattice Q* is not in general equal to that of the direct 
lattice Q as we will see in the next section by explicit computation. 



F. Connection with other formulations 

The definitions we introduced in this section treat the direct lattice and the dual lattice 
separately, using a discrete Hodge operator defined exclusively on either the direct lattice or 
on the dual lattice. This was possible due to the reduction of the problem to a tilted square 
lattice, that can be seen as approximately self-dual. 

In the liter at ure 20 ~ 25 many examples can be found in which the Hodge operator on 
d-dimensional lattices is defined as an operator taking p-cochains on the lattice to (d-p)- 
cochains on the dual lattice. For two dimensional problems, where the Hodge operator is 
closely related with the concept of holomorphy, such operators have been extensively studied 
by Mercat 26 ' 27 , who defines the Hodge operator in terms of the double lattice, the direct sum 
of the direct lattice G and the dual lattice G*. 

In the double lattice formulation the number of holomorphic differentials is double of 
the continuum analogue and we lose either the interpretation of the period matrix as the 
periods of a holomorphic basis or its expected dimension. A 2g x 2g period matrix can be 
defined 27 using a combination of (5.26) and (5.27) together with a procedure that yields a 
g x g period matrix. 
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For a square, pn = 0, version of the lattices G that we study in this paper the f2 and 
Q* period matrices we introduced are closely related with two sub-matrices of the 2g x 2g 
period matrix of the double lattice formulation. The same is not true for general triangular 
lattices G. Since we require the number of holomorphic differentials to be g, and use an 
exact mapping for discrete harmonic differentials from G to G a , our discussion refers to the 
dual of the effective square lattice G n and not to the dual of the original lattice G. 

VI. NUMERICAL EVALUATION OF THE PERIOD MATRIX 

The formalism developed on the previous sections can be illustrated by the numerical 
evaluation from first principles of the period matrices in (3.5). 

We consider lattices with sizes (M 1 , M 2 , K, N 1 , N 2 ) = L(m 1: m 2 , k, n 1: n 2 ) for fixed aspect 
ratios (mi, m 2 , k, ni, n 2 ) and increasing L. The numerically evaluated period matrices for 
both the direct and the dual lattice, fl and Q*, were found to be symmetric with positive 
definite imaginary part for all lattice sizes and coupling constants. A typical example of the 
dependence of the entries of Q and Q* on the lattice size M is shown in Fig. 10. 

The two period matrices were found to have equal real parts for all couplings and aspect 
ratios tested while the imaginary part is typically different but converges towards the same 
value in the large Af limit. In the example of Fig. 10 we see that the imaginary parts of Q 
and Q* converge slower to the continuum limit than their real part, while the average of the 
imaginary parts converges as fast as the real part. 

The period matrix characterizing a given lattice is not uniquely defined but depends on 
the choice of the basis of the first homology group, {aj,6j} used for its evaluation. Period 
matrices evaluated with different bases are related by modular transformations. Consider 
the basis {oj, hi] shown in Fig. 1, we will call this basis HI and H2 to the basis obtained 
from this one by exchanging the labeling of the Oj and hi cycles. 

In table II we give examples of f2 and evaluated at finite A~ and the corresponding 
Af — > oo extrapolation, using the two basis HI and H2. In table III we plot the ratios of 
theta functions Ri for these period matrices, they are related by a permutation of the index 
labeling the half-integer characteristic. From this table we see that modular invariance 
is broken at finite size by Q and Q* but it is satisfied by the average (Q + Q*)/2 with a 
precision of 1% or better. In the large Af limit both f2 and f2* satisfy modular invariance 
approximately while the average period matrix satisfies it with a precision of 10~ 6 . 

Similar averages have been considered in the literature 27 but the fact that the average of 
the two period matrices is modular invariant at finite size and provides a faster convergence 
to the continuum limit is an unpredicted numerical finding. 

The comparison between the averaged period matrix and the period matrices obtained by 
fitting the ratios of determinants in (3.5), is done in table I. There is a remarkable agreement 
between the two sets of values. The theta function ratios corresponding to the two period 
matrices agree with a precision of 10~ 3 or better. A typical example of the convergence of 
ratios of determinants and ratios of theta functions for Q and Q* is shown in Fig. 3. 
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VII. CONCLUSIONS 



In this paper we studied a shape and topology dependent finite size correction to the 
Ising model partition function on genus two triangular lattices. 

We found that, in a similar way to the toroidal case 13,14 , the determinants of the Kaste- 
leyn adjacency matrices converge at criticality to a common bulk term times theta functions 
of half-integer characteristic. This result (3.5) together with the Kasteleyn expression for 
the higher genus partition function (1.2) suggests that the shape dependent correction (1.4) 
is universal for the Ising model on orientable two dimensional lattices. Different topologies 
correspond to sums over different genus theta functions while different lattice aspect ratios 
and coupling constants change only the period matrix Q. The analysis done in this paper for 
genus two can be readily generalized for arbitrary genus since both the Kasteleyn formalism 
and the period matrix evaluation are defined independently of the lattice genus. 

The period matrix characterizing the lattice continuum limit can be extracted from the 
finite size correction (1.4). This provides a unique testing ground for discrete formulations 
of holomorphy and Riemann surface theory. We used such a formulation to reproduce 
the dependence of the period matrix on the lattice aspect ratios and coupling constants. 
The present discussion and the recent interest 26 ' 27 ' 33 in discrete holomorphy may provide 
a starting point for a deeper mathematical understanding of the relation between discrete 
models and certain aspect of conformal field theory. 

We also emphasize that the observed dependence of the determinants of adjacency ma- 
trices on the Kasteleyn orientations is precisely the dependence of the determinant of the 
Dirac operator on the Riemann surface spin structures 10 ' 11 of the corresponding conformal 
field theory. Such a one to one correspondence between discrete lattice determinants and 
functional determinants is reminiscent of the Ray-Singer theorem 34 ' 35 relating determinants 
of the analytical Laplacians on a surface with the determinants of combinatorial Laplacians 
on a triangulation of that surface. We believe that further progress can be made along these 
directions. 
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FIGURES 




FIG. 1. The genus two triangular lattice. The exterior boundary has torus-like boundary 
conditions, the higher genus topology is due to the additional handle obtained by identifying the 
edges numbered with the c and d numbers. The a>i,bi cycles form a basis of the first homology 
group and should be seen as drawn over the lattice edges. All lattice faces are triangular except 
for two hexagons whose edges are marked with the dotted line. 
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FIG. 2. The Ising decorated lattice in the ^4(0000) clockwise odd orientation. The orientation 
of the hexagonal decoration is shown on the inset. The remaining clockwise odd orientations can 
be obtained from this one by inverting the orientations of the edges crossed by a given choice of 
the di, bi cycles. 
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FIG. 3. Comparison between the nine non-vanishing ratios of determinants of adjacency matri- 
ces (full line and star symbols) and ratios of theta functions (dotted and dashed lines and triangle 
symbols) with period matrices evaluated by the procedure of section V. Ratios are shown as func- 
tion of the number of lattice vertices M, for a lattice with aspect ratio (mi,m2,k,ni,n2)=(4,2,2,2,4) 
and couplings [wh,w v , Wd]=[0. 3685,0.1586,0. 2821]. The period matrix evaluations with triangle up 
and triangle down symbols correspond respectively to evaluations done on the direct and on the 
dual lattice. The dashed line is the theta ratio for the average of the two period matrices. 
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FIG. 4. The six vanishing ratios of adjacency matrices determinants for the lattice of Fig. 3. 
These ratios correspond to theta functions of odd half-integer characteristic. 



23 




24 



A 

A, 



-l 



A 

A 



-l 



-l 



-l 



A 

B 



A 

B 



















































I 




















-1 
















-1 





































































FIG. 6. The closed differentials A^jB^. The differentials are zero on all edges except the ones 
crossed by the respective a& loop. These edges are shown in bold on the figure, together with 
the value that the differential takes at that edge. For convenience we use a 62 loop different but 
homologically equivalent to the one shown in Fig. 1. 
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(a) 



(b) 



FIG. 7. The aspect ratio r and tilt angle 9. a) the combinatorial description of the lattice, b) 
the geometrical description of lattice from the parametrization of sinh 2Ki . 
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FIG. 8. The dual lattice Gq is shown, in dashed lines, together with the direct lattice G n , on 
full lines. Each circle is a dual lattice vertex. The two up triangles correspond to a single dual 
lattice vertex, equally for the two down triangles. The numbers labeling the edges that enter the 
octagonal faces show the order in which they join the corresponding vertex. 
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FIG. 9. The dual lattice Gq schematically represented as Gn with two vertices identified on 
each octagon. On the right the local arrangement of the dual edges on the octagonal faces is shown. 
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FIG. 10. The period matrices O and f2* and their average as a function of the number of lattice 
vertices J\f, for the lattice B considered on Figs. 3 and 4. 
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TABLES 





A: [0.426,0.291] (2,2,2,2,2) 


B: [0.369,0.159] (4,2,2,2,4) 


C: [0.082,0.291] (4,2,2,4,2) 


i 


R\ 


flf(fit) 


? 4 (eval) 




f 


ftfit) 


9i(eval) 




flf(fit) 


0?(eval) 


1 


0.15448 


0.15445 


0.15368 


0.01344 





01344 


0.01341 


0.14675 


0.14677 


0.14512 


2 


2.710~ 6 








1.510 -6 










8.6IO- 7 








3 


2.210" 6 








2.010" 7 










4.410~ 6 








4 


2.210" 6 








1.610 -6 










1.010 -6 








5 


5.710~ 7 








6.710" 8 










2.410" 7 








6 


0.87017 


0.87020 


0.87184 


0.35920 





35921 


0.36085 


0.89475 


0.89477 


0.89528 


7 


0.00048 


0.00048 


0.00045 


0.00005 





00005 


0.00005 


0.00117 


0.00110 


0.00111 


8 


0.76230 


0.76227 


0.76487 


0.37222 





37222 


0.37388 


1.00000 


1.00000 


1.00000 


9 


2.710" 6 








5.210" 7 










3.210" 6 








1U 


0.15448 


0.15445 


0.15367 


1.00000 


1 


00000 


1.00000 


0.36001 


0.36003 


0.36012 


1 1 

1_ J. 


0.18531 


0.18534 


0.18381 


0.03254 





03254 


0.03222 


0.10945 


0.10944 


0.10852 


1 9 


0.18531 


0.18534 


0.18380 


0.96246 





96245 


0.96271 


0.32775 


0.32773 


0.32835 


13 


5.710~ 7 








6.210~ 7 










1.510 -6 








14 


0.87017 


87020 


87183 


0.81775 





81775 


81676 


69481 


69483 


69485 


15 


0.03959 


0.03956 


0.03877 


0.03799 





03799 


0.03764 


0.04184 


0.04181 


0.04119 


16 


1.00000 


1.00000 


1.00000 


0.78572 





78572 


0.78496 


0.62338 


0.62336 


0.62345 




fin 


^12 




fin 




«12 


^22 


fin 


^12 


^22 


O^(fit) 


0.2536 


-0.1837 


0.3675 


0.3697 




0.2754 


1.2506 


1.1846 


-0.9445 


1.9637 


n*(fit) 


1.4042 


-1.1978 


2.3956 


0.6491 




0.5826 


2.4293 


1.1483 


-1.0984 


2.5132 


ft r (eval) 


0.2549 


-0.1855 


0.3711 


0.3706 




0.2752 


1.2439 


1.1847 


-0.9429 


1.9583 


n*(eval) 


1.4069 


-1.2031 


2.4063 


0.6495 




0.5841 


2.4343 


1.1480 


-1.1009 


2.5217 



TABLE I. Comparison of ratios of determinants of adjacency matrices with ratios of theta 
functions for the critical Ising model in three different lattices: A, B and C. Each lattice is charac- 
terized by the couplings and aspect ratio [u^,u>,,](mi, 771,2 ,k, 774,712). The L — > 00 ratios are compared 
with theta function ratios for a fitted period matrix J) (fit) and the period matrix f2(eval) evaluated 
by the procedure of section V (first homology group basis HI). 
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first homology group basis: HI 



N = 1008 


^11 ^12 ^21 ^22 


^11 ^12 ^21 ^22 


n 
n* 


0.370493 -0.275048 -0.275048 1.24336 
0.370493 -0.275048 -0.275048 1.24336 


0.643674 -0.566694 -0.566694 2.38225 
0.655324 -0.601642 -0.601642 2.48709 


J\f = CO 


^11 ^12 ^21 ^22 


^11 ^12 ^21 ^22 


n 

n* 


0.370551 -0.275225 -0.275225 1.24389 
0.370551 -0.275225 -0.275225 1.24389 


0.648086 -0.579928 -0.579928 2.42195 
0.650837 -0.588183 -0.588183 2.44671 


first homology group 


r>o qi c • H9 


M = 1008 


^11 ^12 ^21 ^22 


^11 ^12 ^21 ^22 


n 

n* 


-0.871520 -0.216570 -0.216570 -0.220063 
-0.871520 -0.216570 -0.216570 -0.220063 


1.43105 0.334292 0.334292 0.396879 
1.43385 0.328678 0.328678 0.408106 


Af = oo 


^11 ^12 ^21 ^22 


^11 ^12 ^21 ^22 


n 

n* 


-0.871523 -0.216564 -0.216564 -0.220075 
-0.871523 -0.216564 -0.216564 -0.220075 


1.43209 0.33221 0.33221 0.401042 
1.43275 0.330885 0.330885 0.403693 



TABLE II. The period matrix for the lattice B, evaluated for two choices (HI and H2) for the 
basis of the first homology group. Results are shown for both the direct lattice period matrix Q 
and the dual lattice period matrix f2* at finite size N = 1008 and in the thermodynamic limit 
extrapolation. 
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M = 


1008 






N -- 


= 00 






n 


(n + n*)/2 


n 


(n + n*)/2 


HI 


H2 


HI 


H2 


HI 


H2 


HI 


H2 


HI 


H2 


of 


of 


0.013898 


0.013413 


0.013402 


0.013407 


0.013525 


0.013411 


0.013409 


0.013409 


oi 


of, 


0.352897 


0.372173 


0.360822 


0.360606 


0.358967 


0.363532 


0.360852 


0.360852 


Oj 


y io 


0.000068 


0.000031 


0.000049 


0.000050 


0.000053 


0.000044 


r\ r\r\r\r\ ac\ 

0.000049 


0.000049 


ol 




0.366408 


0.385027 


0.373849 


0.373640 


0.372113 


0.376527 


0.373884 


0.373884 


Ofo 


e* 


1.000000 


1.000000 


1.000000 


1.000000 


1.000000 


1.000000 


1.000000 


1.000000 


Oil 


oi 


0.034803 


0.030549 


0.032203 


0.032251 


0.032809 


0.031807 


0.032216 


0.032216 


0\2 


oi 


0.959525 


0.964607 


0.962725 


0.962670 


0.961978 


0.963155 


0.962706 


0.962707 




0l5 


0.811513 


0.812151 


0.816897 


0.816887 


0.815469 


0.815670 


0.816764 


0.816763 




0\A 


0.040773 


0.036186 


0.037616 


0.037667 


0.038358 


0.037285 


0.037639 


0.037639 


o\% 


Ofe 


0.777163 


0.782368 


0.785101 


0.785040 


0.783083 


0.784349 


0.784960 


0.784959 



TABLE III. Comparison of the theta function ratios for the period matrices of table II. The 
ratios corresponding to the period matrix evaluated with the H2 basis, are permuted to make 
modular invariance explicit. 
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